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We report Cu and La nuclear magnetic resonance (NMR) measurements in the title compound that 
reveal an inhomogeneous glassy behavior of the spin dynamics. A low temperature peak in the La 
spin lattice relaxation rate and the "wipeout" of Cu intensity both arise from these slow electronic 
spin fluctuations that reveal a distribution of activation energies. Inhomogeneous slowing of spin 
fluctuations appears to be a general feature of doped lanthanum cuprate. 
PACS Numbers: 76.60.-k, 74.72.Bk, 75.30. Ds, 75.40. Gb 
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Lanthanum cuprate, the prototypical single layer high 
temperature superconductor, has been extensively stud- 
ied for several years to understand the origin of its un- 
usual normal state behavior as well as the mechanism 
for superconductivity. Rare earth co-doped lanthanum 
cuprate has received attention recently because elastic 
neutron scattering experiments have revealed ordering 
of doped holes into charged stripes that constitute anti- 
phase domain walls producing incommensurate antiferro- 
magnetic (AF) order in the intervening undoped domains 
. Charge stripe order is likely intimately related to the 
high temperature superconductivity |^-||. Isostructural 
lanthanum nickelate demonstrates clear stripe order 
and it has been shown there that both the charge or- 
der and the magnetic order are glassy 1^,0. R is also 
known that the magnetic order associated with charge 
ordering in lanthanum cuprate is glassy 1^,^, but this 
situation is more difficult because the charge superlat- 
tice peaks are very hard to observe, presumably because 
the stripes tend to be dynamic. As a consequence little 
detail is known about the glassy behavior. Hunt et al. 
have observed suppression of the Cu NQR signal inten- 
sity ( "wipeout" ) with decreasing temperature that they 
attribute to charge stripe order [lO| . 

NMR provides information complementary to neutron 
scattering because the nuclei are sensitive to the local 
magnetic field and the dynamic behavior of the electronic 
system without requiring spatial correlations. Chou et 
al., first proposed that the very strong peak in the 
^^^La nuclear spin relaxation rate i^Qji-i ^j^g vicin- 
ity of T = 10 K is associated with spin freezing [pd]| . 
Kataev et al. have observed slow spin fluctuations in 
Lai.8-a;Euo.2Sr2;Cu04 at low T Our single crystal 

studies of i^^Tj""'^ have shown that the peak is due to 
continuously slowing electronic spin fluctuations; in par- 
ticular the characteristic fluctuation frequency dis- 
plays an activated temperature dependence [l3| . Fur- 
thermore, these data demonstrate a distribution P{Ea) 
of activation energies Ea centered at Ea/ksT ^ 50 K 
and with a width comparable this center value indicat- 
ing strongly inhomogeneous magnetic properties p3| . To 



understand if this inhomogeneity arises from disorder due 
to, e.g., substitutional dopants, we have applied this anal- 
ysis to several lanthanum cuprate systems exhibiting AF 
order at low temperatures to allow us to explore the ef- 
fect of varying the density and character of the disor- 
der: in-plane doping by Li substitution for Cu, variation 
of doping density in LTT phase Lai.8_a;Euo.2Sr2:Cu04: 
0.01 < a; < 0.15. Remarkably, we find that the character 
of the inhomogeneity, that is, the distribution of activa- 
tion energies is essentially unchanged in all these cases 
and very similar to lightly doped La2_3;Sr2;Cu04 PH , 
suggesting that this inhomogeneity is intrinsic rather 
than arising from impurity disorder. 

These inhomogeneous slow spin fiuctuations also en- 
hance the Cu spin relaxation rate ^^Tj^^and directly ex- 
plain the wipeout of Cu intensity reported by Hunt et al. 
p0| . Because Cu nuclear moments experience a hyperfine 
coupling Am to these fiuctuations that is two orders of 
magnitude larger than for ^^^La and ^ oc A^f, ^Sy-i 
becomes so fast as to relax the signal before it can be 
observed. The distribution of electron spin fiuctuation 
frequencies means that different Cu spins will move out 
of the observation window of the spectrometer at differ- 
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ent temperatures; using P{Ea) obtained from 
quantitatively explain the loss of Cu intensity. 

Co-doping with J — Q Eu rather than Nd is advan- 
tageous: neither the magnetism in the Cu02 plane nor 
the nuclear magnetism suffer the effects of the large Nd 
magnetic moment. The crystal used in this study was 
grown using the traveling solvent floating zone method 
under oxygen pressure of 3 bar Diffraction data as 
well as La NMR indicate a sharp LTT structural phase 
transition at 135 ± 2 K. The observation of incommen- 
surate magnetic peaks near 30 K in Eu co-doped com- 
pounds by elastic neutron scattering |l^ reveals static 
magnetic stripe order. From dc magnetization measure- 
ments, no superconducting transition is observed down 
to 4.2K. /iSR studies of this sample reveal static spin 
order below 25K jl^ . 
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FIG. 1. T dependencies of ^^'^t^^ in lanthanum cuprate 
doped by various routes (Lai.986Sro.oi4Cu04 (o) from pH). 
Solid lines are fits as described in the text using the pa- 
rameters shown in Table 1. Inset: C^^TiT)"^ versus T in 
Lai.65Euo.2Sro.i5Cu04 (solid squares) and Lai.85Sro.i5Cu04 
(open squares; from |r^). 

sition. The spin lattice relaxation rates were measured 
by monitoring the recovery of the magnetization after an 
inversion pulse, and the Cu echo decay was observed by 
integrating the spin echo. The Cu spin lattice relaxation 
and echo decay data were measured in a field of 84.2 kC, 
and the La spin lattice relaxation data were obtained in a 
field of 59.7 kC The intensities of the Cu and La signals 
were obtained from field swept spectra at constant fre- 
quency; the areas A under the spectra were then adjusted 
for echo decay and spin lattice relaxation effects. 

isgji-i jg jjjgai fQj- probing slow electron spin fluctua- 
tions because La is located outside of the Cu02 planes, 
and so only weakly coupled to the electronic system. For 
T > 50 K, except for a range of 35 K around T^t, the 
magnetization recovery is well fit by the standard ex- 
pression for magnetic relaxation of the central transition 
of a spin 7/2 nucleus with a single component Ti 
Between 125 K and 160 K and for T < 50 K, the magne- 
tization recovery was fit with the stretched exponential 
form: M{t) ^ Ma[l - 2 exp(-^i/Ti)], where Mq is the 
equilibrium magnetization. This expression represents 
the magnetization recovery for a distribution of relax- 
ation rates, with a peak at I/6T1. The data for ^^'^T^^ 
are shown in Fig. |l|. The peak at 135 K reflects the LTT 
transition, where the spin lattice relaxation is dominated 
by quadrupolar components |p^ . 

Recently, Suh and coworkers [|l3| demonstrated that 
the strong low temperature peak in '^^^T^^ is accurately 
described by the Bloembergen, Purcell and Pound 

TABLE I. Parameters describing distribution of activation 
energies: P{Ea) = N exp[-{Ea- Eof /2^^]. Fits (see Fig. 
were performed with fixed Too =0.03 ps. 
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mechanism (BPP) Ig^l introduced to explain nuclear spin 
relaxation that results when the characteristic electron 
spin fluctuation frequency (t~^) decreases continuously 
with decreasing temperature: Tc = Tocexp{Ea/kBT). 
The relaxation rate is given by the spectral density of 
fluctuations (typically a Lorentzian) evaluated at the 
Larmor frequency ojl pl|: 



l/T,=-f'hiTj{l + u;iT^) 



(1) 



where ho is the local field fluctuating at the nuclear site, 
and 7 is the gyromagnetic ratio. Hence, the peak oc- 
curs at the temperature where the continuously slow- 
ing characteristic frequency matches the measurement 
frequency: = lol. Since the peak temperature 

is inherently probe-frequency dependent the dynamics 
are best described by i?a- A crucial difference from 
the standard BPP model is that a distribution of Ea 
is required to describe the relaxation data. Roughly 
speaking, the high temperature side of the peak in 
lagji-i determines the center of the distribution, while 
the slow decrease of ^^^Tj^^on low temperature side can 
only be explained if some fraction of the sample expe- 
riences much smaller values of Ea- We choose a Gaus- 
sian as a convenient distribution of activation energies: 
PiEa) = J\fexp[-{Ea - Eo)'^/2A^], where Af is a nor- 
malization factor; we now consider i39y-i func- 
tion of position describing relaxation in some region of 
space at a particular time. This gives rise to a distribu- 
tion of spin lattice relaxation rates P{'^^^T^'^)d^^^T^^ = 



P{Ea)dEa- The measured ^''^T^ is then given by: 
i39y-i ^ J- I39y-ip(i39y-i)^i39y-i^ ^he solid line 

in Fig. |l| is a fit to this expression. 

To explore the dependence of this inhomogeneous dis- 
tribution of magnetic properties on impurity disorder 
we show i39jT|~i results from several systems in Fig. H: 
the same Eu-doped cuprate with an order of magnitude 
lower hole doping: Lai.785Euo,2Sro.oi5Cu04 Jl^ ; lightly 
in-plane doped La2Cu.98Li.02O4 [P2| ; and the Chou re- 
sults ||ll[]: Lai.982Sro.oi8Cu04. The evident similarity of 
the low temperature relaxation data is confirmed by the 
similar distributions of Ea (see Table 1). This indicates 
that impurity disorder is not the crucial element for this 
inhomogeneity, rather it appears intrinsic. 

We observe the same wipeout of Cu signal intensity 
(see Fig. ^ discussed extensively by Hunt et al. [|lO|. A 
is proportional to N/T, where N is the number of nuclei, 
and T is the temperature, therefore N AT is propor- 
tional to the number of nuclei which give rise to the NMR 
signal. In Fig. |2| we show N{T) for Cu and La in various 
fields and orientations in Lai.65Euo.2Sro.i5Cu04. About 
96% of the signal from the Cu nuclei vanishes between 90 
and 4 K, whereas none of the La signal is lost. The Curie- 
Weiss broadening of the spectra is typical for planar Cu 
in cuprates, however the spectra lose intensity uniformly, 
and develop no anomalous features. An apparent loss of 
signal due to a shift of intensity to another region of the 
spectrum can be ruled out. A field dependent study 
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FIG. 2. NMR measurements in Lai.65Euo.2Sro.i5Cu04 
showing the wipeout of Cu signal. The solid symbols rep- 
resent Cu data: (squares) 77 MHz, Ho-Lc; (circles) 86 MHz, 
Ho-i-c; (triangles) 95 MHz, Ho||c. The open diamonds repre- 
sent La data for Jfo||c at 35.9 MHz. The solid line is a plot of 
Io{T) as described in the text with k=1.5. Inset: Ho\\c'''^Cu 
echo size versus the pulse spacing r is plotted for a series of 
temperatures. The fits (solid lines) are described in the text. 

revealed no evidence for a component with a different 
magnetic shift. Even a large change in the local NQR fre- 
quency would be easily detected in the frequency of the 
NMR central transition, where the quadrupolar shift is 
second order (i.e., a 10 MHz change in i^g would shift the 
central transition by ~ 1-2 MHz). However, the temper- 
ature dependence of the Hq _L c quadrupolar shift with 
Hq ~ 80 kG reveals no change in vq with temperature 
down to 4 K. Therefore the evidence suggests that the 
loss of intensity is not due to an inhomogeneity in either 
the local magnetic shift or the electric field gradient. 

Rather, the slow spin fluctuations responsible for the 
strong low temperature peak in ^^^T^^ straightforwardly 
explains the intensity data if we realize that the hyperfine 
field ho of the electronic moments is two orders of mag- 
nitude larger at the Cu site than it is at the La site (see 
Eq. |l|). Thus for nuclei experiencing slow fluctuations, 
Tc > Toos'^ the nuclei relax so fast the spin echo signal 
has decayed before it can be measured; k is a constant 
determined by the recovery time — 15 /isec of the spec- 
trometer following a spin-echo excitation pulse. The echo 
decay of a Cu nucleus experiencing a particular local fluc- 
tuation time Tc is given by: Mq exp[— t/T2_R, — i^/2r|g], 
where T'^ = {(3 + Rf^T{^ is the Redfield term 
and T2G is the Gaussian part of the echo decay as- 
suming static neighbors. Here /3 = 3 and R is the Ti 
anisotropy ratio. Since is distributed, ^^ji-i jg 
well. Intensity is determined by extrapolating the spin 
echo decay data (available for t > 2tr) to t = 0; this 
extrapolation is dominated by the Redfield term, so we 
ignore the Gaussian contribution in the following dis- 
cussion. The full signal from all of the nuclei is then 
given by: M{t) = f Moexp[—t/T2R{Tc)]P{Tc)dTc, where 
P{T,)dTc = P{Ea)dEa. For T > Eo/kB, M{t) decays 



exponentially with the single time constant T2(tc = Too)- 
For T < Eo/kB some regions relax rapidly causing M{t) 
to drop sharply for short t] single exponential behavior 
is then recovered at longer t. One can show that the 
extrapolated intensity loiT) — J^^''^'^ P{Ea)dEa, where 
K = ln(Tcut/'''oo)- For T > Tcut excessive signal decay oc- 
curs for t > 2tr, thus wiping out the signal [we use the 
criterion M{2tr;Tcut) = 1%]- The measured low tem- 
perature anisotropy of i^^Tj"^ |l^ shows ~ hence 
i? ~ i. Using Eq. ^ and assuming A\^f = 100 kOe//iB 
(the wave vector dependence of the very slow spin fluc- 
tuations is not known) gives Tcut = 0.13 ps =^ k=1.5. In 
Fig. H we show a plot of Iq (T) using this value of k and 
the distribution P{Ea) extracted from fits to the La Ti 
data; the agreement of this straightforward model with 
the data is good, especially given the uncertainty in pa- 
rameters. This strongly indicates that Cu wipeout is not 
a measure of the stripe order parameter pO| | . 

The Cu spin lattice relaxation rate, ^'^Tj^^, was mea- 
sured between 4 and 300 K. The inversion recovery data 
were fit to the standard expression for magnetic relax- 
ation of the central transition with a single Ti ||l^ . The 
temperature dependence of (^^TiT)~^ is shown in the in- 
set of Fig. |l|, as well as data from Ohsugi et al. , for 
Lai.85Sro.i5Cu04. The spin lattice relaxation rate of the 
Cu that contribute to the NMR signal at low tempera- 
tures reveals no unusual behavior, and is in fact quite 
similar to the (^^y^j-)-! Lai.85Sro.i5Cu04. Note that 
(63TiT)-i is identical for the two systems for T > 50 K. 
Below this the Lai.85Sro.i5Cu04 becomes superconduct- 
ing at 38 K, whereas the Lai.65Euo.2Sro.i5Cu04 does not. 

The character of the Cu echo decay (Fig. || in- 
set) also provides evidence for Ti inhomogeneity. 
The temperature dependence of the echo decay in 
the Lai.65Euo.2Sro.i5Cu04 is similar to that in other 
cuprates, except for a distinct crossover from Gaus- 
sian to exponential decay around 100 K. In spite of 
the distributed Redfield contribution discussed above 
we would expect Gaussian behavior at large t; this 
is not observed experimentally. This behavior can be 
understood if the measured Cu nuclei are coupled to 
neighboring spins that experience a fast spin lattice re- 
laxation. For like neighbors the echo decay is given 
by M{t) = Moexp[-t^f{t/T^)/2Ti^], where f{x) = 



8x-^[5x/2 + 9e~"=/2 - 2e-^ - 7] l23|. For infinite Ti the 
echo decay is Gaussian; as Ti gets shorter, the echo de- 
cay becomes exponential. The solid lines through the 
data in the inset of Fig. || are fits in which we assume 
y-i ^ S3yi~^cas+®^^iTcxtra> whcrc the cxtra contribution 
is a variable parameter that represents the contribution of 
fast relaxation on the neighbors . Qualitatively, these 
fits indicate that the observed Cu nuclei are coupled to 
neighbors that undergo Ti fluctuations much faster than 
those of the observed nuclei. 

Combined, the La and Cu data reveal an inhonioge- 
neous, glassy freezing of the spin dynamics in the heav- 
ily doped, stripe ordered systems. La relaxation demon- 
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strates a dramatic slowing of spin fluctuations such that 
their characteristic frequency matches ~ 5-50 MHz 
around 10 K, and the data are well modeled by a distribu- 
tion of activation energies. The quantitative explanation 
for the loss of Cu intensity from spin freezing provides 
a direct demonstration that the fluctuations are inho- 
mogeneous: at any given temperature < T < 100 K 
a T-dependent fraction of the Cu spins experience spin 
fluctuations so slow (see Eq. |l|) that they are "invisible" 
while the remaining spins remain visible due to much 
faster spin fluctuations. 

We have seen that quantitatively similar inhomoge- 
neous freezing of spin fluctuations (that are well mod- 
eled by very similar distributions of activation energies) 
occurs in a variety of hole doped lanthanum cuprates con- 
taining impurities of very different character: i) lightly 
Sr-doped lanthanum cuprate containing only a few out- 
of-plane impurities, ii) lightly Li-doped La2Cui_xLix04 
where the impurities are in the Cu02 planes, iii) 
both lightly and heavily Sr-doped Lai.65Euo.2Sro.i5Cu04 
where Eu co-doping adds more out-of-plane impurities 
and induces the LTT structure thought to pin charged 
stripes. The fact that the spin freezing is independent 
of both impurity density and the location of the impu- 
rity (in or out of the plane where the fluctuating mo- 
ments and the charged stripes reside) strongly suggests 
that extrinsic impurities are not an essential element of 
the inhomogeneity, rather it appears to be an intrinsic 
response of the hole doped system to perturbation. Fi- 
nally Hunt et al. have shown a similar loss of Cu in- 
tensity in La2-a;Sr2:Cu04 throughout the doping range 
< X < further emphasizing the generality of 

this phenomenon. We note that loss of NMR signal has 
been observed previously in classical spin glass systems 
p5| and in stripe-ordered lanthanum nickelate [ p6| . 

Two different mechanisms have been explored in the 
literature in order to explain the spin freezing in the anti- 
ferromagnetic state: one involves the "cluster spin glass" 
and freezing of the transverse spin degrees of freedom 
PI ; another the loss at low temperatures of collective 
hole/stripe motion In the presence of frustration a 
Heisenberg system exhibits glassy behavior of the trans- 
verse spin components below the longitudinal ordering 
temperature [ p7| , and Kivelson et al. have pointed out 
the role of disorder in striped systems The ob- 

served very slow spin fluctuations imply large clusters of 
spins (presumably frustrated by disorder) are involved. 
The role of charged stripes in these dynamics isn't clear. 
Stri pe motion is certainly a mechanism for spin dynam- 
ics |22[ , however if the same mechanism is responsible 
for the inhomogeneous relaxation over the wide regime of 
hole density we observe, then the very weak hole density 
(hence stripe density) dependence of the distribution of 
activation energies implies stripe-stripe interactions play 
no role over a decade of stripe density; this seems un- 
likely. Because the stripes constitute anti-phase domain 
walls in the antiferromagnet, stripe disorder leads to se- 
vere spin disorder [BIR]. Hasselmann and Castro-Neto 



have pointed out that for a; ~ 1/8 the frustration can 
arise from point defects and disorder in the stripe topol- 
ogy The slow inhomogeneous dynamics could arise 
from frustrated interactions between such large clusters. 
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